I . INTRODUCTION 


^"Measurements  of  the  ambient  ocean  noise  for  sea  states 
greater  than  zero  show  that  the  noise  is  angular  dependent.  This 
implies  that  the  spatial  gain  of  a beam  forming  sonar  is  dependent 
on  the  orientation  of  the  sonar  array.  There  is,  therefore,  an 
optimum  gain  available  for  a given  sonar  system  operating  in  a 
particular  environment  which  is  orientation  dependent.  This  op- 
timum orientation  can  be  predicted  from  theory  providing  the 
ambient  noise  field  properties  are  known.  However  such  is  not  the 
case.  To  date  there  is  ample  data  in  the  frequency  band  from  0.5 
kHz  to  1.5  kHz  but  very  little  data  above  this  frequency.  In 
addition,  environmental  conditions  such  as  bottom  and  surface  con- 
ditions, velocity  profiles,  inhomogeneities,  etc.  are  not  well 
chronicled  for  the  measured  data.  However  the  latter  situation 
is  not  critical  if  one  were  interested  in  finding  average  values 
for  the  ambient  noise  by  measurements  over  a long  period  of  time 
since  these  differences  average  out.  This  approach  is  very  costly, 
however,  since  a large  quantity  and  variety  of  data  is  required. 

Another  approach  to  the  problem  can^be  taken  which  is 
cheaper  and  can  give  estimates  of  the  ambient  noise  in  agreement 
with  large  data  set  averages.  This  approach  recognizes  that  the 
far  lEield  noise  anisotropy  is  basically  due  to  the  surface  dis- 
turbance and  attempts  to  calculate  the  surface  source  from  data  — 
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as  a function  of  frequency,  sea  state  (or  wind  scale)  and  angle. 
However  the  observed  noise  data  and  its  environmental  conditions 
which  are  used  to  calculate  the  surface  sources  must  be  well  known 
to  insure  that  the  source  function  represents  the  real  world. 

Thus  one  can  predict,  quite  accurately,  the  noise  anisotropy  for 
a particular  environment  if  one  knows  the  source  corresponding  to 
the  surface  disturbance.  This  source  can  be  calculated  from  a 
completely  different  set  of  boundary  conditions  since  the  source 
is  uniquely  related  to  the  sea  state  conditions  only.  TRACOR  has 
taken  the  latter  approach  to  predicting  the  ambient  anisotropic 
noise.  This  document  is  an  offshoot  of  some  of  this  work  done  in 
conjunction  with  systems  performance  studies.  There  is  now  a re- 
quirement to  establish  a library  of  source  data  versus  sea  state 
as  c function  of  frequency  and  angle  so  that  systems  modeling 
problems  for  a large  variety  of  conditions  can  be  solved. 


II.  THE  DIRECTIONAL  AMBIENT  NOISE  MODEL 


The  directional  ambient  noise  field  is  assumed  to  be 
the  superposition  of  two  noise  fields, 

N (SS,  f,  Qr,  «S)  = Nx  (f)  + n2  (SS,  f,  0R,  <*)  (1) 

where  0R  = elevation  angle  of  the  received  noise  in  the  vert- 
ical plane,  0°  along  the  horizontal  axis,  and 
positive  angles  toward  the  sea  bottom. 

i>  - azimuthal  angle  of  the  received  noise,  positive 
angles  being  measured  counterclockwise. 

N^  is  a noise  field  that  is  isotropic  and  which  corresponds  to 
that  which  would  be  measured  in  a sea  state  zero,  perfectly  calm, 
deep  ocean.  Its  magnitude  is  solely  a function  of  frequency. 
^2(®R’^)  assumec*  to  originate  from  a uniform  distribution  of 
noise  generators  which  lie  at  the  sea  surface. 

Its  magnitude  and  shape  will  be  a function  of  frequency, 
f,  depth,  h,  and  ocean  environment.  Measurements  show  that  it  is, 
to  within  about  1 dB,  independent  of  azimuthal  angle  but  this  de- 
pendence can  be  handled  by  the  theory  if  desired.  The  total  di- 
rectional noise  field  for  a particular  sea  state  neglecting  (6 
dependence,  is 


N (0p , f,  h)  *=  N,  + N9  (0p,  f,  h)  . 
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Suppose  that  a receiver  is  placed  at  a depth  h in  the 
deep  ocean  (Fig.  1).  Let  Eq  equal  the  noise  power  per  unit  solid 
angle  being  emitted  from  a point  source  on  the  ocean  surface. 

Let  dAQ  be  the  area  intersected  by  a cone  of  rays  which  subtend 
the  solid  angle  at  unit  reference  distance.  Then  the  energy 
crossing  the  area  dAQ  per  unit  time  and  contained  within  the  cone 
is 


I dA 
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where  I = the  intensity  at  unit  reference  distance  of 

the  surface  element  in  the  direction  0 , <p  . 
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The  assumption  of  ray  acoustics  is  that  the  small  cone 
of  energy  leaving  the  point  source  element  will  diverge  or  con- 
verge as  the  rays  trace  out  their  respective  paths  but  the  total 
power  transmitted  is  constant  for  a constant  number  of  rays.  In 
this  case  the  intensity,  I,  at  the  point  R where  the  receiver  is 
located,  must  satisfy  the  equation 


or 
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where  dAR=  incremental  area  of  the  cone  of  rays  at  the 
receiver,  which  corresponds  to  the  total 
number  of  rays  in  the  solid  angle  described 


by  dAQ. 
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The  spreading  loss  from  a point  S on  the  surface,  to 
a point  R,  at  the  receiver,  is  defined  as  (see  Figure  1) 


so  that 
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Anderson,  Gocht,  and  Sarota,—  have  shown  that  the  spread- 
ing loss  from  a point  S to  the  point  R is  related  to  the  backward 
spreading  loss  from  R to  S along  the  same  path  by 
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where  = velocity  of  sound  at  the  receiver. 

Cg  = velocity  of  sound  at  the  sea  surface. 

The  intensity  of  the  received  noise  may  now  be  written  in  terms 
of  the  received  angle  and  the  spreading  loss  from  the  receiver 
to  the  surface  element  as 
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From  Physics  of  Sound  in  the  Sea—  the  spreading  loss  from  R to  S 


can  be  shown  to  be 
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where  x = horizontal  surface  range 
0R  = angle  at  the  receiver 
0g  = angle  of  ray  at  the  sea  surface 

The  noise  power  propagated  from  a point  source  surface 
element  will  also  suffer  absorption  in  its  journey  to  the  receiver, 
so  that  the  intensity  received  at  angle,  0D,  along  the  direct  path 
(no  surface  or  bottom  bounces)  is 
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where  s = direct  path  length  to  surface  element 
a = attenuation  coefficient 


Equation  (9)  gives  the  intensity  received  via  the 
direct  path,  at  angle  0R,  from  a group  of  point  source  noise 
radiators  at  the  sea  surface.  A more  realistic  assumption,  how- 
ever, is  to  assume  that  the  noise  originating  at  the  sea  surface 
consists  of  a continuous  distribution  of  noise  generators  dis- 
tributed over  the  ocean  surface.  In  accordance  with  this 


viewpoint,  assume  that  each  inf initesmal  area  S of  the  ocean 
surface  generates  noise  power  per  unit  solid  angle  of  (see  Figure  2) 


where 


dEQ  = a G(M)  dS 
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= surface  noise  density 


= normalized  intensity  of  the  noise 
radiated  surface  element  in  the  di- 
rection (®g  > $g) • 


then  the  noise  intensity  received  is  given  by 
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The  surface  element  dS  (see  Figure  2)  can  be  expressed 
in  terms  of  0R  and  <f)  as 

dS  = x dx  d<i> 


But  x is  a function  of  0R  and  the  depth,  h,  so  that  for  constant 
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The  final  expression  for  dIR(0R , <z$)  becomes  in  terms  of  the 
solid  angle  dQ,  since  &&  is  the  same  for  the  source  and  receiver, 


dIR(9R,rf)  = 
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The  above  development  is  for  those  rays  which  arrive  at 
the  receiver  by  direct  paths.  In  addition  to  noise  arriving  over 
the  direct  path  from  a surface  element  at  horizontal  range  x^, 
noise  may  arrive  at  the  angle  0R  from  surface  elements  located  at 
ranges  x.2>  x3  • • • xi  with  corresponding  path  lengths  of 
s2>  s3>  • • • sj*  (See  Figure  3).  If  such  rays  strike  the  sea 
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surface  or  bottom  they  will  experience  corresponding  energy  losses. 
Accordingly,  for  rays  arriving  from  the  surface  which  have  struck 
both  the  bottom  and  surface,  the  total  noise  received  from  all 
such  elements  is 
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(13) 


Y = bottom  loss  coefficient 
0 = surface  loss  coefficient 


For  those  rays  which  suffer  surface  loss  but  do  not  strike  the 
bottom,  the  noise  intensity  is  given  by  equation  (13)  with  the 
bottom  loss  coefficient  set  equal  to  one.  For  rays  which  have 
no  direct  path  to  the  surface  without  first  striking  the  bottom, 
the  noise  intensity  is  given  by 
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The  total  noise  intensity  per  unit  solid  angle  is  ob- 
tained by  summing  the  sea  state  zero  noise  per  unit  solid  angle 
and  the  appropriate  expression  from  equation  (13)  or  (14) 

N(0R,rf)  = nx(0r)  + N2(9r)  (15) 


The  total  noise  received  from  all  directions  is  ob- 
tained by  integrating  N(9R,$)  over  all  solid  angles  at  the  receiver. 


III.  COMPUTER  PROGRAM 


A computer  program  for  calculating  equation  (15),  the 
noise  per  unit  solid  angle,  and  equation  (16),  the  total  noise 
intensity,  has  been  written  and  perfected  for  the  Univac  1108 
digital  computer.  Additional  features  of  the  program  are  the 
ability  to  solve  for  a G(9),  the  product  of  the  surface  noise 
density  and  the  noise  source  radiator  pattern,  from  experimentally 
measured  data  and  to  give  the  directional  noise  field  over  a fre- 
quency band.  Major  points  of  the  computer  model  are  as  follows. 

For  any  specified  depth,  frequency,  sea  state,  bottom 
loss,  velocity  profile,  and  a G(0g),  the  directional  noise  field 
in  the  vertical  plane  is  calculated  in  increments  of  1°  from  -89° 
(directly  overhead)  to  +89°  (toward  the  sea  bottom) . The  results 
of  this  noise  field  computation  are  then  used  to  calculate  the 
total  noise  intensity  by  performing  a numeric  integration  over 
all  solid  angles  at  the  receiver.  If  experimentally  measured  data 
of  N(9)  is  read  into  the  program,  along  with  pertinent  environ- 
mental and  velocity  profile  data,  the  program  will  solve  for 
0 G(9)  and  give  the  results  in  increments  of  1°.  If  the  direc- 
tional noise  over  a frequency  band  is  desired,  rather  than  single 
frequency  data,  the  directional  noise  field  is  calculated  for  fre- 
quency increments  of  specified  length,  and  the  resultant  noise 
field  is  obtained  by  integrating  over  the  bandwidth. 


IV.  ANALYSIS  OF  EXISTING  DIRECTIONAL  NOISE  FIELD  DATA 


The  computer  model,  based  upon  the  results  of  Section 
II,  has  been  used  to  analyze  the  rather  fragmentary  experimental 
directional  ambient  noise  data  which  now  exists.  The  results  of 
this  analysis  to  be  discussed  below,  give  a better  physical  under- 
standing of  the  directional  ambient  noise  and  give  support  to  the 
general  validity  of  the  model. 

A.  THE  SURFACE  SOURCE  FUNCTIONS 

If  the  surface  noise  source  density  function,  a G(0),  is 
known,  then  equation  (13)  and  (14)  allow  the . calculation  of  the 
directional  noise  field  for  any  specified  velocity  profile,  depth, 
bottom  loss,  surface  loss,  frequency,  and  sea  state  conditions. 
Becken^{  Urick— { Forster—^  and  others  (Ref.  1,  2,  4,  6,  9)  have 
postulated  normalized  G(0)  radiators  of  the  form  cosn9,  where  n is 
an  integer  or  fractional  power,  but  thus  far  no  measured  data  has 
been  used  to  obtain  G (0 ) . If  measured  directional  noise  field 
data,  and  the  associated  physical  conditions  under  which  it  was 
obtained,  are  known,  then  it  becomes  possible  to  solve  for  the 
source  density  function,  a G(0),  from  equations  (13)  and  (14). 
Using  measured  directional  ambient  noise  field  data  taken  by 
Axelrod,  Schoomer  and  Von  Winkle^  empirical  expressions  for  G(6) 
have  been  obtained  from  the  computer  model. 
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The  measured  directional  noise  data  of  Von  Winkle  was 
obtained  with  the  Trident  Array,  which  sits  on  the  ocean  bottom 
at  a depth  of  14,400  feet  about  25  miles  south  of  Bermuda.  Data 
is  available  for  low  frequencies  from  112  to  1411  Hz  and  as  a 
function  of  Beaufort  wind  scale.  The  data  are  averaged  over  a 
period  of  time  from  Winter  thru  Summer,  and  specific  velocity 
profiles  for  a measured  set  of  data  are  not  available.  Average 
velocity  profiles  for  an  area  13  miles  from  the  Trident  site  have 
been  supplied  by  Von  Winkle*  however,  and  appropriate  velocity 
profiles  from  this  data  have  been  used  for  obtaining  the  G(0) 
functions.  Bottom  loss  as  a function  of  grazing  angle  and  fre- 
quency were  taken  from  Marsh^.  Surface  loss  was  assumed  to  be 

unity  for  the  low  frequencies  under  consideration,  and  the  atten- 
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uation  as  given  by  Thorp—'  was  used.  Under  these  conditions  the 
normalized  G(0)  noise  source  functions  for  frequencies  of  560, 

891,  1122,  and  1411  Hz  and  Beaufort  sea  state  5 as  obtained  from 
the  Von  Winkle  data  are  given  in  Figures  4,  5,  6 and  7. 

B.  PREDICTED  DIRECTIONAL  AMBIENT  NOISE  FIELD  CHARACTERISTICS 

The  empirically  derived  a G (9 ) functions  as  discussed 
above  have  been  used  to  study  the  properties  of  the  directional 
noise  field  as  a function  of  depth,  velocity  profile,  bottom  loss, 
and  frequency.  Computed  noise  fields,  using  the  empirically 
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derived  a G(0)  functions,  show  that  the  shape  of  the  noise  field 
is  a highly  dependent  function  of  velocity  profile  structure, 
depth  and  bottom  loss.  Typical  results  are  presented  below. 

1 . Depth  Dependence 

Figures  9 through  14  show  a typical  variation  of  the 
noise  field  with  depth,  for  a frequency  of  1411  Hz,  the  velocity  pro- 
file shown  in  Figure  8,  and  Marsh's  bottom  loss.  The  primary  change 
in  the  shape  of  the  noise  field  as  depth  varies  occurs  around  zero 
degrees  elevation  angle.  As  the  depth  at  which  the  noise  field  is 
computed  is  varied,  the  noise  arriving  ac  the  nearly  horizontal 
angles  is  greatly  dependent  upon  the  sound  velocity  at  the  re- 
ceiver relative  to  the  velocity  at  the  ocean  surface.  When  the 
velocity  at  the  surface  is  greater  than  that  at  the  receiver  depth, 
noise  from  surface  will  not  be  able  to  arrive  at  the  receiver  at 
the  nearly  horizontal  angle,  i.e.,  there  exist  shadow  zones  at 
the  surface  which  do  not  contribute  to  the  noise  field.  If  the 
velocity  at  the  receiver  is  equal  to  or  greater  than  the  surface 
velocity  then  there  will  be  no  shadow  zone  areas  at  the  surface 
and  hence  no  null  will  appear  at  the  nearly  horizontal  angles . 
Greatest  variations  in  the  noise  field  shape  as  depth  changes 
occur  in  approximately  the  first  3000  feet  of  depth,  where  the 
velocity  profile  is  changing  most  rapidly. 
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2 . Velocity  Profile  Dependence 


A comparison  of  Figures  9,  16  and  18  along  with  the  ac- 
companying velocity  profiles  of  Figures  8,  15  and  17  gives  an  in- 
dication of  the  velocity  profile  dependence  of  the  noise  field. 

The  depth  of  the  fields  is  near  the  surface  in  each  case.  It 
appears  that  the  characteristic  shape  of  a noise  field  for  a 
specified  profile  is  largely  dependent  upon  the  profile  structure. 
The  sharp  peaks  which  arise  on  both  sides  of  the  null  in  some 
cases  have  a physical  explanation.  These  peaks  arise  because  the 
noise  arriving  per  unit  solid  angle  is  coming  from  a very  large 
surface  area,  due  to  the  near  horizontal  aspect  angles  of  the 
sound  rays  leaving  the  sea  surface  near  the  shadow  zones.  Low  sea 
state  data  taken  by  G.  Fox-'  with  Trident  array  tend  to  confirm 
the  physical  reality  of  such  peaks. 

3 . Bottom  Dependence 

The  lower  half  of  the  directional  ambient  noise  field 
curves  is  strongly  dependent  upon  bottom  loss  as  a comparison  of 
Figures  10  and  19  will  reveal.  Both  curves  are  computed  for  the 
same  sets  of  conditions  except  for  bottom  loss.  In  Figure  10  the 
theoretical  bottom  loss  of  Marsh  was  used.  In  Figure  19  a bottom 
loss  as  published  by  B.  A.  Becken  was  used.  The  bottom  loss 
affects  the  shape  and  magnitude  of  the  lower  half  of  each  curve. 

In  addition  the  magnitude  of  the  upper  half  of  each  cuve  is 
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affected  for  those  rays  which  have  arrived  after  striking  the 
bottom  but  the  effect  is  much  less  pronounced  for  the  upper  half. 


4.  A Comparison  of  Computed  Directional  Noise  Fields 
with  the  Measured  Becken  Data 

The  only  known  measured  data  for  directional  noise  field 
near  the  surface  are  the  four  curves  published  by  B.  A.  Becken  in 
his  thesis—' . The  exact  physical  conditions,  such  as  velocity 
profile,  bottom  loss,  etc.,  under  which  the  Becken  data  were  ob- 
tained are  not  known  and  therefore  exact  reproduction  of  the  Becken 
results  is  not  possible.  By  making  certain  assumptions  as  to  the 
conditions  under  which  the  Becken  data  were  taken  it  is  possible 
to  show  a general  qualitative  agreement  between  the  measured  data 
of  Becken  and  the  calculated  noise  using  ct  G(0)  calculated  from 
Von  Winkle's  observed  noise. 

The  Becken  data  are  not  single  frequency  data  but  are 
taken  over  a 750-1500  Hz  band.  The  curves  of  the  restored  or 
"true"  field  that  are  published  were  obtained  from  the  measured 
data  by  attempting  to  remove  the  array  effects  of  the  system  used 
for  measurement.  The  array  used  by  Becken  for  the  noise  measure- 
ments had  a rather  "fat"  main  lobe  beamwidth  and  measurements  were 
taken  only  in  5°  steps.  Any  attempted  restoration  of  the  noise 


field  measured  in  such  a manner  will  omit  much  of  the  fine  detail 
of  the  noise  field  structure^/. 


Because  of  the  above  considerations  it  was  felt  that  the 


easiest  way  to  check  the  measured  results  of  Becken  with  the  com- 
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puted  noise  field  was  to  compare  the  resultant  curves  of  the  noise 
that  would  be  measured  at  the  output  of  a Becken  type  array  for 
both  Becken' s restored  data  and  TRACOR's  computed  noise  field. 

To  do  this  the  Becken  array  as  published  in  his  thesis  was  simu- 
lated on  the  1108  Univac  computer.  Theoretical  noise  field,  for 
frequencies  of  560,  891,  1122,  and  1411  Hz,  were  then  computed  for 
conditions  of  a Pacific  Ocean  velocity  profile  assumed  to  be 
similar  to  that  under  which  the  Becken  data  were  measured.  The 
results  were  then  integrated  over  a frequency  band  from  750  to 
1500  Hz  to  produce  the  curve  shown  in  Figure  20.  Using  this  curve 
to  represent  the  noise  field  in  the  water  and  the  comparable 
"corrected"  data  were  integrated  over  Becken' s array  beam  pattern 
to  obtain  the  noise  out  of  the  array  as  a function  of  elevation 
angle.  The  results  are  presented  in  Figure  21. 
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V.  AREAS  FOR  FUTURE  INVESTIGATION 


There  remain  many  improvements  to  be  made  on  the  de- 
veloped noise  model  and  further  investigation  needs  to  be  under- 
taken to  determine  the  characteristics  of  the  directional  noise 
field.  Some  of  these  are  as  follows: 


1.  For  higher  frequencies  a surface  loss  should  be  in- 
corporated into  the  model. 

2.  Improved  accuracy  could  be  obtained  by  incorporating 
continuous  velocity  profile  rather  than  straight  line  approxima- 
tions for  calculating  the  path  length. 

3.  The  contribution  of  shadow  zone  areas  to  the  noise 
field  should  be  investigated.  At  the  present,  sharp  disconti- 
nuities occur  in  the  noise  field  for  the  angles  at  which  the  shadow 
zone  occurs.  Results  published  in  Physics  of  Sound  in  the  Sea 
indicate  that  the  sound  falls  off  as  a sharp  exponential  from  the 
shadow  zone  and  that  there  will  be  some  contribution  to  the  noise 
intensity  from  such  zones.  Presently  shadow  zone  contributions 

are  considered  negligible. 

4.  Solutions  for  the  a G(9  ) surface  source  density 

s 

noise  functions  need  to  be  studied  further.  Because  of  velocity 
profile  structure  it  is  difficult  to  obtain  accurate  solutions  for 


a G ( 9 s)  for  angles  near  zero  degrees.  The  accuracy  for  these 
angles  becomes  important  in  predicting  the  intensity  of  the  noise 
field  for  angles  arriving  nearly  horizontal  at  the  receiver. 

5.  Most  important  of  all,  more  experimentally  measured 
data,  along  with  the  accompanying  physical  conditions  under  which 
it  was  measured,  needs  to  be  obtained.  From  data  measured  under 
the  proper  conditions,  it  will  be  possible  to  obtain  more  reliable 
a G ( 0 ) functions  as  a function  of  frequency  and  sea  state.  Once 
these  have  been  obtained,  it  then  becomes  possible  to  accurately 
model  the  directional  ambient  noise  in  a manner  analogous  to  that 


described  above. 
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